Giant vortices, vortex rings and reentrant behavior in type-1.5 superconductors 
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We predict that in a bulk type-1.5 superconductor the competing magnetic responses of the two 
components of the order parameter can result in a vortex interaction that generates group-stabilized 
giant vortices and unusual vortex rings in the absence of any extrinsic pinning or confinement 
mechanism. We also find within the Ginzburg-Landau theory a rich phase diagram with successions 
of behaviors like type-1— >• type- 1.5— >type-2— >type-1.5 as temperature decreases. 
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Topological defects have been intensively investigated 
in a wide variety of physical systems in e.g. condensed 
matter (ferromagnets, superfluids, liquid crystals, etc.) 
[1|, optics 0, particle physics or astrophysics [3J, which 
has given rise to fruitful exchanges of ideas between the 
different fields. In a superconductor the quantized line 
defect is a vortex (also called flux line) which carries a 
localized integer multiple of the magnetic flux quantum 
$o = h/2e. Close to the transition temperature a unique 
quantity, the Ginzburg-Landau (GL) parameter k, de- 
termines superconductivity response to a magnetic field 
[4j . It is defined as the ratio of the magnetic penetration 
depth A to the coherence length £ of the order parame- 
ter (OP) and it can be calculated from the microscopic 
parameters of the material within the BCS theory. In 
a type-1 superconductor, k < 1/V%, the interaction be- 
tween vortices is purely attractive which results in their 
fusion into macroscopic normal domains in the interme- 
diate state. On the contrary in a type-2 superconductor, 
K > l/v2) the interaction is purely repulsive so the vor- 
tices are stable and form a lattice of co-directed lines 
in the mixed state. At the critical value the GL the- 
ory (which is mathematically equivalent to the abelian 
Higgs model in particle physics) predicts that the vortex 
interaction is exactly zero [5|, |6j. However experiments 
with k f=a l/\/2 observed complex patterns consisting of 
mixed-state domains with vortex-free Meissner domains 
[J, [Z|. This stems from a non-monotonic inter- vortex 
potential which is repulsive at short separation and at- 
tractive at a long range. Such interaction profile can be 
explained theoretically by taking into account different 
mechanisms such as low-temperature corrections 4 , |8| or 
fluctuations and anisotropy in the vortex lattice 8 1 . 

Recently Babaev and Speight [l0( predicted that in 
a superconductor with two uncoupled components the 
inter-vortex potential can be non-monotonic when one 
component of the OP is intrinsicly of type 1 while the 
other is of type 2. This idea has been recently supported 
by observations [ll| of highly inhomogeneous vortex dis- 
tributions in clean samples of two-band superconductor 
MgB2 [12j. The estimates of MgB2 parameters [llj indi- 
cate that this behavior, coined type-1.5 superconductiv- 



ity, only occurs in high quality samples which explains 
why the material was classified before as type-2. Includ- 
ing the finite inter-band Josephson coupling (which is 
moderate in MgFJ2), GL calculations [llj yielded non- 
monotonic vortex interaction potentials, also when one 
of the bands is inactive. A detailed study of this inter- 
band 'proximity effect' [13] has shown that a type-1.5 
potential is then associated with dissimilar variations of 
the OP components at finite distances from the vortex 
core, despite their long-range behaviors being identical. 
Type-1.5 superconductivity could be more commonly ob- 
served than the secondary effects showing up in single- 
component superconductors with k = l/\/2- Stimulated 
by the growing number of discovered multi-band super- 
conductors, e.g. iron- based superconductors [1J, LL5|, we 
investigate in this article novel features of the type-1.5 
regime which distinguish it from single-component su- 
perconductivity at the critical k. We find, in particular, 
group-stabilized giant vortices and vortex rings composed 
of single/multi-quantum vortices as thermodynamically 
stable phases. 

The two-band GL theory. — In the context of two-band 
superconductivity with s-wave symmetryof the OP we 
consider the GL free-energy functional [16j 

Tgl =?x+Ti- 7(*i*2+*i*a) + 7T^ x A ) 2 (1) 
where the contributions from each band 
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are supplemented by a Josephson-type coupling and the 
magnetic energy. The Josephson interaction describes 
Cooper pair tunneling between the two bands. It is the 
simplest of the terms mixing the two components \IJ n of 
the OP and in MgB 2 , for example, the other ones van- 
ish in the clean limit [lg, |l7|. Although a\ and ai can 
change sign at different temperatures, any finite Joseph- 
son coupling 7 imposes that both ^ n vanish at the same 
critical temperature T c . Below T c in the homogeneous 
Meissner state A = and |$ n | = u n , which arc solutions 



of u n (a n + (3 n u 2 n ) - 
then 



7U3_ n = 0. The energy density is 
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Superconductivity with two OP components is richer 
than with a single one because not only the ratio u 2 /u\ 
changes with temperature but also the two components 
can have different magnetic responses and spatial varia- 
tions. Let us introduce the band parameters K n and their 
average R defined by 



K n = $ m ?lV //V(27r) 3/2 
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where the condensation-energy fraction 
p n = (3 n u^/2\e M \- 



(4) 
(5) 



(6) 



If the two components were separated not in band space 
but in real space as in a Josephson junction, k\ and k 2 
would be the actual GL parameters of the two supercon- 
ductors in contact. Note that R varies with temperature 
contrary to K n . To discuss the different types of vor- 
tex interaction it is convenient to rescale the distance 
x = X^/2x' , the vector potential A = ($o/"\/27rA)a and 
\fr„ = u n ip n . Here the magnetic penetration depth is 
defined by 



A 2 = $o«/(47r)^|e M | 



3 / 2 U»J 1 /2 



(7) 



The normalized functional T = (Fgl — £m)Aa is then 

Z7T 



+ ri|D^i| 2 + r 2 \D^ 2 \ 2 + -(V x a) 2 } 



(8) 



where the line energy e\ — (<f'o/47rA) 2 , D — (— iV + a), 
the normalized Josephson coupling factor 

q = -iuiu 2 /\e M \ (9) 

and the supcrfluid-density fraction 

r„ = M 2 J m~ 1 /( u i m i~ 1 + u l m 2 1 )- ( 10 ) 

So in contrast to single-component superconductivity 
where the shape of the vortex interaction depends only 
on k ||, the shape depends here on four parameters: 
{R,pi, q, ri} (since p\ +p 2 = 1 and r\ + r 2 = 1) or equiv- 
alently {R, k\, k 2 , q}. However the Josephson term locks 
ipi with tp 2 when q diverges to infinity at the critical tem- 
perature (since u n vanishes) and, as illustrated below, the 
single-component behavior with -01 = 02 and k — R is 
recovered in this limit if R 7^ l/v2- 

Numerical evaluation of the vortex interaction. — The 
vortex-pair potential is numerically calculated by varia- 
tional minimization of the GL energy F — jFd 2 x with 




FIG. 1: (color online) Potential of the vortex interaction 
U vv (d) = F(d) — F(oo) as a function of vortex separation d for 
different values of average GL parameter R (band parameters 
ki = 0.5, K2 = 5) and normalized Josephson coupling q. 



the constraint that the pair separation d is held fixed [18j . 
The Ansatz [6] for F is built upon isolated-vortex solu- 
tions with winding number L = 1 and 2. For the latter 
the radial distributions ipn{ r ) an( i a ( r ) are found by solv- 
ing the finite difference version of the GL equations on a 
linear grid of 3000 points with a relaxation method. The 
boundary conditions are \%pi\ = \^ 2 \ — a# = at r = 
and |^i I = I ^2 1 = 1, a$ = —L/\/2r at large distance 
r - I 00A. 

Fig. Q] shows the interaction potential U vv (d) = F(d) — 
F(oo) of two vortices when one OP component is of type- 
I («i = 0.5) and the other one of type-2 (n 2 = 5) while 
the average R is 0.6, I /^/2 and 1. For R — 0.6 {R = 1) the 
energy of a two-quanta vortex is smaller (larger) than the 
energy of two isolated vortices like in a single-component 
type-1 (-2) superconductor. But we note that this is not 
always the case for other choices of parameters. Besides 
for small q, the inter-vortex potential is repulsive at short 
separation and attractive at long range even though R is 
not in the vicinity of \j\/2. As q increases the potential 
is first of type 1.5 and then it recovers a type-1 (-2) shape 
while the position <i m i n of its minimum goes to (+00). 
In contrast, the potential for R = l/v2 is always of type 
1.5 whatever the strength of the coupling. 

Distribution of vortices. — To illustrate the variety of 
magnetic responses in a type-1. 5 superconductor Fig. [2] 
shows the vortex patterns resulting from different interac- 
tion potentials U vv ranging from the more type-1-like (i) 
to the more type-2-like (hi). They are obtained by molec- 
ular dynamics simulations of a system of over-damped 
vortices (see Ref. ll| for details). Starting from several 
random distributions the dynamics of 200 vortices was 
simulated in a square of size 200A x 200A with periodic 
boundary conditions during 100000 time steps [19j . 

The highly inhomogeneous vortex distributions are 
characteristically composed of areas filled with vortices 
surrounded by or surrounding vortex-free regions, de- 
pending on the ratio of the minimal-potential separation 
d m in to the average distance d vv between vortices (com- 



pare e.g. (b-iii) with (b-iv) in Fig. [5]). When d min > d v 
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FIG. 2: (color online) (a) Different potential of vortex inter- 
action in a type-1.5 superconductor. The inset is a magnifica- 
tion of curve it), (b) The vortex distributions resulting from 
the interaction potentials i to iv. L<I>o-vortex means a vor- 
tex with L quanta of flux. The dashed circle represents the 
positions of U vv minimum around one vortex. 



the whole surface is covered by a vortex lattice even 
though the interaction is of type-1.5. With potentials 
i and U a multi-quantum vortex can be stabilized by 
the surrounding vortices at distances ~ d m in: the lat- 
ter generate a local potential well which overcomes the 
repulsion between the merged vortices forming the giant 
vortex. However the maximum of U vv at zero separa- 
tion prevents all the vortices from collapsing into one 
unique macroscopic vortex, i.e. a normal region like in 
the intermediate state. Giant vortices can be observed 
in mesoscopic samples [20J or trapped by pinning sites 
2l| . But the ability to stabilize them in the bulk with- 
out any extrinsic confinement mechanism is one unique 
property of type-1.5 superconductivity. Yet we found 
an even more striking feature: rings of vortices can ap- 
pear, as shown in Fig. [31 when the OP components vary 
on length scales which are greatly different (i.e. when 
n\ <C k-2 and q -C 1). In this situation the local max- 
imum of the inter-vortex potential at zero separation is 
sharper than the minimum (see curves i and v in Fig. [2] 
and [3j which is detrimental to the stability of a giant 
vortex. The favored configuration is then an annular dis- 
tribution of radius < d m i n , with the presences of one cen- 
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FIG. 3: (color online) (a) Vortex-interaction poten- 
tials for {k, Ki,K2,q} = {0.65,0.5,25,0.4} (curve i) and 
{0.83,0.5,100,0.08} (curve t>). (b) Different vortex rings ob- 
tained with potential v and resulting from local variations of 
the vortex density (same conventions as in Fig. [2]). 



tral vortex and of multi-quantum vortices if the vortex 
density is high. 

Limit of large Josephson coupling. — The robustness 
of type-1.5 superconductivity at R« l/v2 towards the 
large coupling q can be analytically demonstrated. When 
q is large the difference (ipi —^2) is small so the GL func- 
tional ([8]) can be minimized with a perturbative expan- 
sion in powers of q^ 1 by writing ip±.2 = fa ± ipi/2q. 
At zeroth order in <? _1 , (fo and a are then solutions 
of the Bogomolnyi-Jacobs-Rebbi (BJR) equations when 
k sa 1/V2 [a |(| |8[ while the first-order correction ipi = 
(P-R/2R 2 )(l-\<po\ 2 )vo with 



P = Pi - P2 and R = ri - r 2 . 



(11) 



At zeroth order the free energy does not vary with the 
distance between the vortices since the energy of BJR 
solutions does not depend on the positions of the vor- 
tices. Its dependence on the separation is given at first 
order of correction by the variation of Jd 2 x[(R 2 — 1/2 + 
2s)|<po| 4 — slvol 6 ] which can be estimated in the lim- 
its of short and large separation as in Ref. [8|. Here 
s = (2k 2 P-R) 2 /4qk 2 . We find H3 that for large Joseph- 
son coupling q the vortex interaction is 

of type 2 when C = ^Zr/S^ > °- 346 ' 
of type 1.5 when -0.4 < C < 0.346, 

of type 1 when C < -0.4. (12) 

P and R quantify the disparities between the two bands. 
When they are small or when q is large, as in the vicin- 
ity of T c , the ratio C grows out of the bounds defining 
type 1.5 and superconductivity is either of type 1 or of 
type 2. However in the case R = 1/V2, C is exactly 
zero which means that the pair potential is always non- 
monotonic independent of the strength of the Josephson 
coupling. 

Temperature dependence of the interaction. — Since 
the average R and the normalized Josephson coupling 
q depend on temperature the vortex-vortex interaction 
can also change with it. As an illustration, the tempera- 
ture evolution of these quantities and of the ratio C are 
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FIG. 4: (color online) (a) Temperature evolution of nor- 
malized Josephson coupling q, average GL parameter R and 
vortex-interaction parameter C for {ki,K2J = (i) {0.5,1.5}, 
(ii) {0.6,2}, (iii) {0.6,2.65}, (iv) {0.65,3} and (v) {0.7,4}. 
(b) Temperature evolution of the separation d m i n at U vv - 
minimum and of the magnetic penetration depth A for the 
parameter set (iii) in (a). 



plotted in Fig. Ufa) for different values of {k±, K2} which 
yield type- 1.5 superconductivity away from T c . We used 
parameters which are in the range of experimental esti- 
mates for MgB 2 [Uliilli!]. We ad °P ted the model of 
temperature dependence for a n from Ref. lg, [23 1 with 
«i becoming negative below T± — 0.64T C and «2 below 
T-2 = 0.91T C . At the critical temperature q diverges to in- 
finity so the vortex-vortex interaction is either of type 1 
or type 2 just below T c . For the majority of param- 
eter sets there is only one change of behavior type-l(- 
2)— >type-1.5 as temperature decreases. But as shown by 
the example (iii) in Fig. [4^a) R can pass by the value 
1/V2 in the vicinity of T c while q is large. In this case 
C can go out of the bounds (|12[) twice which means an 
alternation of four regimes is possible. It was previously 
found [ll|, [l3[ that a type- 1.5 phase can survive above 
min(Xi,T2) but a second interval above max(Ti,T2) was 
unexpected. This four-phase scenario is deduced from in- 
equalities (fl"2"|) that are valid only when q is large. Thus 
U vv for the set of parameters (iii) is numerically eval- 
uated to check the estimation. We find that the phase 
diagram is indeed composed by an unusual reentrant suc- 
cession of four regimes: type-1— >■ type-1. 5— >type-2— ^type- 
1.5. The temperature dependence of the position d m i n 
and of the penetration depth A are plotted in Fig. HJb). 
When the superconductor enters a type-1 phase the po- 
sition of the non-monotonic-potential minimum goes to 
zero while it diverges to infinity when entering a type-2 
phase. 

In conclusion, we predict novel properties for type-1. 5 
superconductivity which strikingly distinguish it from the 
vortex states of conventional type-1 and type-2 super- 
conductors. By numerical GL calculations of the vortex 



interaction potential and simulations of the resulting vor- 
tex dynamics, we found vortex-stabilized giant vortices 
which survive in bunches and unusual vortex rings in bulk 
multiband superconductors. The multi-component char- 
acter produces as well a complex phase diagram where 
type- 1.5 behavior can re-enter by alternating with type- 
1 and type-2 interludes in a succession of three regime 
transitions as temperature is changed. 
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EF/05/005 (INPAC) from the University of Leuven. 
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